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The low-energy theory of d-wave quasiparticles coupled to 
fluctuating vortex loops that describes the loss of phase coher- 
ence in a two dimensional d-wave superconductor at T = is 
derived from first principles. The theory has the form of 2+1 
dimensional quantum electrodynamics (QEDz), and is pro- 
posed as an effective description of the T = superconductor- 
insulator transition and of the pseudogap phase in under- 
doped cuprates. The coupling constant ("charge") in this 
theory is proportional to the dual order parameter of the XY 
model, which is assumed to be describing fluctuations of the 
phase of the superconducting order parameter. Finiteness of 
the charge is then tantamount to the appearance of infinitely 
large vortex loops, i. e. to the loss of phase coherence in 
the system. The principal result is that the destruction of su- 
perconducting phase coherence in the d-wave superconductors 
typically, and immediately, leads to the appearance of antifer- 
romagnetism. This transition can be understood in terms of 
the spontaneous breaking of an approximate "chiral" SU C (2) 
symmetry, which may be discerned at low enough energies 
in the standard d-wave superconductor. The mechanism of 
this spontaneous symmetry breaking is formally analogous 
to the dynamical mass generation in the QED3, with the 
"mass" here being proportional to staggered magnetization. 
Other phases with broken chiral symmetry include the trans- 
lationally invariant "d+ip" and "d+is" insulators, and the 
one- dimensional charge-density and spin-density waves, which 
are all insulating descendants of the d-wave superconductor. 
All the insulating states have the neutral spin-1/2 excitations 
that one can identify in the superconductor confined by the 
logarithmic potential. Electron repulsion is in this formal- 
ism represented by a particular quartic perturbation to the 
QED3 action, which breaks the chiral symmetry and selects 
the antiferromagnet as the preferred broken symmetry state. 
I formulate the mean-field theory of the antiferromagnetic in- 
stability in presence of a short-range repulsive interaction, and 
find the staggered magnetization to be significantly enhanced 
deeper inside the insulating state. The theory offers an expla- 
nation for the rounded d-wave-like dispersion seen in ARPES 
experiments on the insulating CaiCuOiCli (F. Ronning et. 
al, Science 282, 2067 (1998).) Relations to other theoretical 
approaches to the high-T c problem are discussed. 

PACS: 74.20.Mn, 74.25.Jb, 74.40.+k 



I. INTRODUCTION 

Soon after the original discovery, it became well ap- 
preciated that the high temperature (high-T c ) supercon- 
ductors are all quasi two-dimensional insulating antifer- 
romagnets that become superconducting with the intro- 



duction of holes. The nature of the relationship between 
antifcrromagnetism and high temperature superconduc- 
tivity has been the central issue in the field. Follow- 
ing the time honored strategy of understanding first the 
non-superconducting state, most of the approaches to 
the high-T c problem focused on finding the mechanism 
by which doping an antiferromagnet would produce a su- 
perconductor H . The essential difficulty in pursuing this 
strategy seems to be that the Mott insulator is itself a 
non-trivial strongly correlated state, harder to describe 
in simple terms than the metallic Fermi liquid, which 
played its role in the BCS theory of the low-temperature 
superconductivity 0] . The situation becomes only worse 
away from half-filling, where the ground state of even the 
simplest models becomes more ambiguous. Experimen- 
tally, the cuprates seem to loose their antiferromagnetic 
ordering with doping before they become superconduct- 
ing, and many candidates for the intermediate "pseudo- 
gap phase" have been discussed in literature. The na- 
ture of the non-superconducting state that is supposed 
to be unstable to superconductivity with doping is at 
this point, however, far from clear, and may prove to 
be non-universal. Arguably, the physics of underdoped 
regime may be the main mystery of high temperature 
superconductivity. 

In a remarkable contrast to the uncertainties inher- 
ent to the insulating phase, the superconducting phase 
of most high-r c materials is well established to have the 
d-wave symmetry of the order parameter ||, [Q, typ- 
ically with well-defined, long-lived quasiparticle excita- 
tions ||, ||. This simplicity suggests that an inverted 
approach to the high-T c problem may be more natural 
[]7j: if there exists a d-wave state in the phase diagram, 
which other states can in principle be inferred from it? 
The purpose of this paper is to establish the theoretical 
framework for answering this question, answer it, and 
show how this may help explain some salient features of 
the cuprates phase diagram and the angle resolved photo 
emission spectroscopy (ARPES) experiments in the in- 
sulating state ||, |)[. 

Loosely speaking, there are two ways to destroy a su- 
perconducting state: 1) by driving the amplitude of the 
order parameter to zero, which is what is well described 
by the weak-coupling BCS theory at finite temperature 
Q, for example. For a d-wave superconductors this pro- 
cess presumably is relevant at large dopings, where weak- 
coupling treatments of the Hubbard and related models 
can be trusted, and disorder should eventually force T c 
to vanish |u|. 2) Even if the amplitude of the order pa- 
rameter is large and finite, superconductivity will be lost 
with the destruction of phase order jn|, p2[ . There is 
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evidence that this is what actually occurs in underdoped 
cuprates, where the superconducting transition tempera- 
ture (T c ) is much lower that the pseudogap temperature 
T* . Since underdoped cuprates are strongly two dimen- 
sional, at finite temperatures the loss of phase order may 
be expected to proceed via the Kosterlitz-Thouless tran- 
sition, and indeed, there are distinct experimental sig- 
natures of the fluctuating vortices above T c [^4[ . 
The following question then naturally arises: What is 
the nature of the T = phase that derives from a 
two-dimensional d-wave superconductor when the phase- 
coherence is lost, but the order parameter amplitude is 
still finite? The central thesis of this work is that phase 
incoherent d-wave superconductor (dSC) is nothing but 
the insulating (typically incommensurate) spin-density- 
wave (SDW), i. e. weak antiferromagnet. Short account 
of this result appeared earlier in Q . 

I show that the minimal continuum theory of the low- 
energy quasiparticle excitations near the four nodes of 
the d-wave order parameter coupled to fluctuating vortex 
loops at T = is provided by the 2 + 1 dimensional 
quantum electrodynamics (QED 3 ): 

S = J d 2 fdr[^ a ,(d, + *a,)*< + ^pyp( V x C 1 ) 

where v = 0, (imaginary time) 1, 2 (space), and the sum 
over repeated indices is assumed. The four- component 
Dirac fields $j j = 1,2 represent the sharp, electrically 
neutral spin- 1/2 excitations one can define in the super- 
conducting state (and hence may call "spinons"), which 
are minimally coupled to a massless gauge-field a. The 
gauge-field derives from the fluctuating topological de- 
fects (vortex loops) in the phase of the superconducting 
order parameter, which have been integrated out in deriv- 
ing the theory (1). Complex number ($) is proportional 
to the the disorder (dual order) parameter Jl6| , and rep- 
resents the state of vortex loops: (<&) ^ signals the ap- 
pearance of infinitely large loops in the system and the 
loss of phase coherence, which is the T = analog of the 
Kosterlitz-Thouless transition |l7|] . In the superconduct- 
ing state, on the other hand, all loops are of finite size, 
($) = 0, and the gauge field, in the simplest approxima- 
tion, may be considered effectively decoupled from the 
spinons: quasiparticle excitations are then sharp, since 
all the short-range interactions that have not been explic- 
itly written in the Eq. (1), if weak enough, are strongly 
irrelevant. When ($) ^ the situation becomes radi- 
cally different, as the gauge- field mediates a long-range 
interaction between spinons. In reality the theory is also 
strongly anisotropic, but for simplicity this possibly im- 
portant feature has been neglected in writing the Eq. 
(1). The QED% has also been recently considered by 
Franz and Tesanovic |lj] as an effective description of 
the pseudogap state. They argued that the presence of 
the massless gauge field may explain the broad features 
seen in ARPES measurements in the normal state 
Jl9| . Here I show that at T — as soon as ($) becomes 
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FIG. 1. A schematic phase diagram of cuprate supercon- 
ductors in terms of the low-energy chiral symmetries SU C (2) 
(full) and U c (l) (broken). Besides the chiral symmetries, 
d-wave superconductor (dSC) also has the full spin rotational 
symmetry, and the spin-density wave (SDW) has the super- 
conducting £7(1) and the spin rotational symmetry around 
one of the axis. Near and left to the underdoped transi- 
tion point the system is proposed to be an extremely weak 
SDW, which becomes reinforced by the electron repulsion, 
and which continuously evolves into a stronger antiferromag- 
net near half-filling. 

finite there is a dynamical generation of the mass term 
~ m^i^i in (1), which can be identified as the stag- 
gered potential felt by the original electrons, i. e. with 
the SDW order parameter. Quantum fluctuating dSC is 
thus at T = inherently unstable towards SDW ordering 
once the phase coherence is lost. 

The dSC — > SDW quantum phase transition is an ex- 
ample of spontaneous breaking of a continuous global 
symmetry in the Eq. (1), which for a lack of better name 
I will call "chiral" throughout the paper. Chiral symme- 
try breaking is a well studied field-theoretic phenomenon, 
believed to be inextricably linked to confinement in the 
QED3 pC[ |. Massless QED3 for single species of Dirac 
fermions has the continuous £7(2) = £7(1) X SU C (2) sym- 
metry, with the generators /, 73,75, and 735 = 17375, re- 
spectively. In the action in the Eq. (1), the £7(1) factor 
represents the residual spin rotational symmetry left by 
the choice of representation, as will be explained in de- 
tail later. It is the additional SU C (2) symmetry per Dirac 
component in the QED3 that will be of central interest 
here. The fermion mass term mM>i^> breaks the SU C (2) 
for each Dirac field to £7 C (1), and the two broken gener- 
ators rotate between different insulating states. Chiral 
SU C (2) arises as an approximate symmetry of the dSC 
only at low-energies, and will be manifestly broken, for 
example, by higher order derivatives omitted in the Eq. 
(1). It should not be confused with the spin rotational 
symmetry which is, of course, also, and exactly, present 
in the dSC. Higher order derivatives and the electron in- 
teraction terms reduce the SU C (2) to its £7 C (1) subgroup, 
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which is related to the spatial translations of the original 
electrons. The identification of the approximate chiral 
symmetry in the dSC is essential for establishing the con- 
nection between the antiferromagnetic and the supercon- 
ducting phases advocated in this paper, and represents 
one of my central results. The idealized cuprates phase 
diagram may be understood in terms of the chiral sym- 
metries of different states as depicted in Fig. 1. 

Assuming the scale for the SDW transition T$dw{x) in 
an anisotropic quasi-two dimensional high-temperature 
superconductor to be set by the magnitude of the stag- 
gered magnetization at T = |HJ , the present work sug- 
gests that near and left of the superconductor-insulator 
transition one should expect it to be considerably lower 
than the superconducting T c (x) near and right of the crit- 
ical point: Tsdw{%u — 5) << T c (x u + 5), where x u is the 
critical doping for the dSC-SDW transition, and S << 1 
(see Fig. 1). This is because the generalized QED3 with 
TV fermion species has a critical point at iV = N c « 3, 
above which there is no dynamical mass generation [ p0[ . 
The QED 3 in (1) has N = 2 components, which together 
with some numerical factors gives very weak SDW order 
near the superconducting phase. The pseudogap phase 
in cuprates at T — is therefore proposed here to be 
actually an extremely fragile SDW, likely to be easily 
destroyed by disorder, for example. As half-filling is ap- 
proached and the vortex loop condensate ($) increases, 
the repulsion between electrons also becomes important. 
Short-range repulsion is represented in the QED3 by a 
particular quartic term, which if weak is irrelevant in the 
superconducting state, but which also manifestly breaks 
the chiral symmetry of the low-energy theory. I show 
that the effect of such a term is first to break the de- 
generacy among states with broken chiral symmetry in 
favor of the SDW, and then to dramatically increase the 
SDW order parameter farther from the dSC. The pic- 
ture implied by the QED3 is qualitatively in accord with 
the generic phase diagram for the underdoped cuprates, 
where the antiferromagnetic transition near half-filling 
raises to ~ 300K, but is typically unobservably low very 
near the superconducting state. 

Neutral spinons, which are well defined quasiparticles 
in the superconducting state, in the insulator become 
broad excitations with the lifetime proportional to the 
antiferromagnetic order parameter. At T = and at 
large distances they become confined by a logarithmic 
potential provided by the gauge-field in presence of the 
chiral symmetry breaking. Due to the weakness of the 
SDW order very near the superconducting transition, 
however, spinon confinement is effective only at very 
large distances, or equivalently, at very low tempera- 
tures. The weak SDW phase therefore appears effectively 
deconfined at intermediate length scales. The finite-T 
pseudogap phase has the gapless spinons strongly scat- 
tered by the massless gauge field, in qualitative agree- 
ment with the broad spectral features of the electrons 
seen in ARPES @. Near half-filling the SDW order in- 
creases and the bound state of spinons rapidly shrinks, 



leaving only magnons in the excitation spectrum. 

Confined nature of the standard antiferromagnet close 
to half-filling, if postulated, by itself already points to 
the QED3 as a viable candidate for the effective the- 
ory of underdoped cuprates. If one views the super- 
conducting state as being spin-charge separated 0, one 
needs a mechanism by which spinons would eventually 
become confined in the antiferromagnetic phase. The 
QED3 provides such a mechanism automatically, since 
the massless gauge-field mediates a long-range logarith- 
mic interaction between the spinons that binds them at 
all energies. Were the gauge-field massive, on the other 
hand, the physics would be equivalent to Z 2 gauge theory, 
and the antiferromagnetic state would be deconfined and 

§'te different from the usual antiferromagnet j22[, |p3| , 
. The very existence of an ordinary antiferromagnet 
at, and presumably near, half filling |25| may therefore be 
taken as evidence in favor of the type of theory presented 
in this paper. 

The physical picture of the antiferromagnetic (SDW) 
insulator as a phase-disordered d-wave superconductor is 
further supported by the ARPES data on the insulat- 
ing Ca 2 Cu0 2 Cl 2 and Sr 2 Cu0 2 Cl 2 §], §. These ex- 
periments show two unexpected features of the insulat- 
ing state: 1) although the ARPES spectral function is 
broad, one can nevertheless identify a remnant of the 
Fermi surface, 2) the dispersion at such an approximate 
"Fermi surface" has a d-wave form, except that it be- 
comes rounded and without the characteristic cusp at low 
energies. The "relativistic" dispersion for broad quasi- 
particle excitations that the QED% implies in the insu- 
lating state, when measured from the lowest energy given 
by the dynamically generated chiral mass, provides a very 
good fit to the data (see Fig. 5). The present theory im- 
plies that the rounding of the dispersion is controlled by 
the size of the sublattice magnetization, and therefore 
should decrease with doping, as one approaches the su- 
perconducting state. It would be desirable to test this 
prediction in future experiments. 

In the body of the paper I develop the above picture in 
detail. In the next section, I derive the Dirac represen- 
tation of the Hamiltonian for low-energy nodal quasipar- 
ticles, and discuss the coupling to quantum fluctuating 
vortex loops in the section III. A derivation of the dy- 
namics of the gauge field starting from the XY model 
on a lattice is presented in the sec. IV. This section is 
somewhat technical and may be skipped at first read- 
ing. Instead, the reader may consult the Appendix B, 
where a simpler derivation for finite temperatures is pre- 
sented. Dynamical breaking of the chiral symmetry and 
the formation of the SDW state is discussed in the sec. 
V. More general discussion of the chiral symmetry and 
the other ordered states on the chiral manifold is pro- 
vided in the sec. VI. The reduction of chiral symmetry 
by the irrelevant terms is discussed in sections VII, and 
the mean-field theory of the antiferromagnetic instabil- 
ity of the QED3 in presence of the electron repulsion is 
solved in the sec. VIII. Confinement of spinons in the in- 
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sulator is discussed in the sec. IX. The discussion of the 
ARPES measurements is given in the sec. X. Summary 
of the main results and the discussion of the relations to 
other theoretical approaches is given in the concluding 
section. I finish with the list of open problems. Techni- 
cal details are presented in five Appendices. 



II. DIRAC THEORY FOR NODAL EXCITATIONS 

I begin by assuming that the superconducting state, 
except from being a d-wave, otherwise exhibits the stan- 
dard BCS phenomenology. In particular, I take that 
the quasiparticles are well-defined, long-lived excitations. 
Generally, the quasiparticle action at T ^ may then be 
taken to be 

S = T ^2 [{iu n - ^)cl(k,uj n )c a (k, u> n ) (2) 

-|A(fe)4(fc, u> n )cl a {-%, -U n ) + h.C. + 0(c% 

where A(fc) has the usual d-wave symmetry, and two spa- 
tial dimensions (2D) are assumed, c and c' are the elec- 
tron operators, a = ± labels the z-projection of electron 
spin, and u) n are the fermionic Matsubara frequencies. 
Units are chosen so that h = c = e = 1. 0(c 4 ) term 
stands for all possible short-range interactions between 
quasiparticles. 

We may represent the quasiparticle Hamiltonian in 
terms of two four- component fields, 



*/(<?*, o>„) = (c\(k,u> n ),C-(-k, -w„), 
c+(k- Qi,u! n ),C-(-k + Qi,-U! n )), 



(3) 



where Qi — 2Ki is the wave vector that connects the 
nodes within the diagonal pair i = 1, 2, as in Fig. 2. For 
spinor 1, k = Ki+q, with \q\ <C \Ki\, and analogously for 
the second pair. The construction of the four-component 
field is not unique. The choice in the Eq. (3) differs from 
the one made in the ref. for example. I postpone the 
discussion of the alternative construction used there for 
the Appendix D. Using the construction in the Eq. (3), 
and by observing that £g = and A £ = -Ag^., 

for k « Ki, and then by linearizing the spectrum as £g = 
Vfq x + 0(q 2 ) and A^ = v&q y + 0(q 2 ), one arrives at the 
low-energy action 

= J d 2 r dr^[d T + M lVf d x + M 2 v A d y ]*' 1 (4) 

+ (1 -> 2, x <-> y) + 0(9*'^*', *' 4 ), 

with (3 = 1/T. The continuous Dirac field ^(r,r) is 
defined as 



d 2 q 
tt^) 2 
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FIG. 2. The wavevectors Ki, i = 1,2, and q. The dashed 
line stands for the putative Fermi surface. The SDW ordering 
wave vectors are Qi = 2Ki. 

with the integral over momenta performed over \q\ < A < 
T*. The 4x4 matrices in the Eq. (4) are Mi = 10-3 ® er 3 , 
and M% — —io~3®o~i. a are the usual Pauli matrices, and 
the coordinate system has been rotated as in Fig. 2. 

To cast the theory in Dirac form we may invoke the 
matrix 70 = o\ (£> I, where / is the 2x2 unit matrix. 
Then -f 2 = / ® /, and M, = 707^, with 71 

= 2<L 



<J2 <X> CT 3 , 

and 72 = -a 2 ®ai. {7^,7^} = 2b~ vll , v,/j, = 0, 1, 2, so the 
7-matrices indeed satisfy the Clifford algebra. The quasi- 
particle action (2) at low energies becomes equivalent to 
the field theory 



d 2 f 



dr^odr + jiVfdx + jiVAdy}^ (6) 
(1 2, x ^ y) + 0(d* W, *' 4 ), 



where — ^//to- Weak quartic interactions, as long 
as they are short-ranged, are irrelevant by simple power 
counting. This simply reflects the severe phase-space re- 
strictions for scattering of the nodal quasiparticles. I will 
therefore omit them temporarily, together with the sec- 
ond order derivative terms, ato return to their effects in 
the section VII. 

The reader would be correct to note that there is a 
considerable freedom in selecting the form of the matrix 
70. In fact, any 4x4 matrix that anticommutes with 
M\ and M 2 and squares to unit matrix would yield an 
equally valid Dirac representation. It is shown later that 
this freedom will correspond to different "directions" in 
the space of ordered states with broken chiral symmetry. 
The specific choice for 70 made here will be analogous to 
choosing a direction in real space along which to search 
for a finite magnetization, for example, in the more fa- 
miliar magnetic phase transitions. 



(5) 



III. COUPLING TO TOPOLOGICAL DEFECTS 

The goal in this section will be to find the most eco- 
nomical form of the coupling between nodal excitations 
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in the dSC and the fluctuations of the phase of the su- 
perconducting order parameter. Working assumption is 
that the amplitude fluctuations are frozen well below the 
pseudogap temperature T* , so it is only the phase degree 
of freedom that remains active at low energies. With this 
in mind I write 

v A v A (r,r) = |u A | e «(<Mr»+Mr») ; (7) 

where <f> r represents the regular (" spin- wave" ) part of the 
order parameter phase, and (f> 8 is the singular contribu- 
tion due to topological defects. At T = these would 
be the vortex loops O], or the more familiar vortices 
and antivortices at T ^ 0. At this point it is tempting 
to transform both spin-up and spin-down fermionic op- 
erators by absorbing a half of the total superconducting 
phase into each. In presence of topological defects, how- 
ever, this would lead to multivalued fermionic fields and 
would not be a local change of variables in the partition 
function. This problem may be circumvented by allowing 
only vortices of double vorticity 0, for example, which 
then leads to the Z 2 gauge theory representation of the 
problem, and a possibility of spin-charge separation in 
the pseudogap regime Q. It is the single vortices, how- 
ever, that first become relevant at the T ^ Kosterlitz- 
Thouless transition S|, and they should be included in 
the description of the T = transition as well. I will 
therefore utilize the idea of Franz and Tesanovic [p7| , 
p8| who suggested dividing a given vortex configuration 
into two groups A and B, and transforming the electron 
operators with spin up and spin down differently. We 
write 

M?,T) = ^f± + <t> sA {r,T), (8) 

and similarly for B. 4> s a is the piece of the singular part of 
the phase that comes from the defects grouped in A. One 
may then make a local change of variables by introducing 
a new Dirac field \& as 

*(f,r) = U^t^'^t) (9) 

where U = diag{e~^, e^ fl , e -4 ^, e^ B }. Since any 
given vortex defect is either in group A or B, and there- 
fore associated either with up, or with down spin by 
the transformation (9), circling around it with the trans- 
formed fcrmion would yield either 2n or zero of the ac- 
cumulated phase change. Components of the new field \& 
are therefore single- valued functions. 

The gauge-transformed action for the Dirac field is 
then 

S[W] -> S[f,a,v\ = (10) 
J d 2 r J dr*i[7o(9 r + ia ) + jiVf(d x + ia x ) 
+72|i>a|(<3j/ + ia y )]^i + (1 —> 2,x <-> y) + iv^J^, 
with a v = d u (4>A - 0b)/2, V v = d v {(j>A + $b)I% and 



Since the vector J„ is built out only of the products of the 
creation and the annihilation operators with same spin, 
it also represents the physical charge current carried by 
the quasiparticles. On the other hand, since the regular 
part of the phase <p r was in the Eq. (8) divided equally 
between spin up and spin down, the Dirac field '3/ is in- 
variant under a regular gauge transformation. Compo- 
nents of \& therefore create electrically neutral excitations 
with spin- 1/2 §, which may therefore be referred to as 
spinous. 

The action (10) has two rather different gauge sym- 
metries, and it may be worthwhile pausing a little to re- 
flect on them. First, the physical electromagnetic gauge 
field would enter the action (10) by the replacement 
v ii — y v n + A^, and couple to the charge current. Un- 
der a regular gauge transformation — > A^ + d^Xi the 
Volovik's field |29| — > — <9 M x> while the gauge field 
a M and 4" remain the same. The action (10) is therefore 
gauge invariant, in the standard sense. But it also has 
an additional internal gauge symmetry, under the trans- 
formation a M —>■ + d^Xi -l V ~~ * v v ^ ~ * e~ lxi b. This 
reflects the freedom of choice in the Eq. (8); one could 
have equally well chosen the regular part of <\>a to be 
(0r/2) + an d 01 <Pb, (0r /2) — X- O ne deals with this 
gauge freedom, as usual, by eventually introducing the 
gauge-fixing term for that allows one to freely sum 
over all regular internal gauges x- Similarly, the division 
of the singular part of the superconducting phase into 
that which comes from the defects in the group A, and 
the defects in the group B, is equally arbitrary. Just like 
one effectively sums over all regular internal gauges by the 
introduction of the gauge-fixing term, we will sum over 
all singular internal gauges by averaging over all possible 
divisions of defects into two groups. This is explained in 
the next section, and in the Appendix B. As a byproduct, 
the averagings over regular and singular internal gauges 
will insure that up and down spinons are treated equally 
in the QED 3 , in respect of the symmetry of the original 
electronic action (2). 

The crucial observation about the action (10) is that 
the coupling of spinons to phase fluctuations is furnished 
by two U(l) fields which play quite different roles in the 
problem. Total superconducting phase determines the 
Volovik's field v v and couples to the charge current, in 
the same way as the true electromagnetic field would. v u 
will therefore inevitably become massive once the high- 
energy spinons in the Eq. (10) are begun to be inte- 
grated out. Its fluctuations therefore may provide only a 
short-range interaction between spinons. The gauge field 
a„, on the other hand, enters (10) in a gauge- invariant 
way, and therefore is protected from acquiring a mass 
from spinons. Both gauge fields, however, depend on 
the fluctuating positions of the topological defects, and 
acquire their dynamics not only from the spinons, but 
from the defects as well. To determine their dynamics 
one therefore needs to integrate the defect degrees of 
freedom out. If a M would stay massless even after this 
integration is performed, it would mediate a long-range 
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interaction between the nodal excitations, which, unlike 
the short-range quartic terms in the Eq. 6, would not 
be made irrelevant by the phase space restrictions. This, 
however, depends on the precise way a M acquires its dy- 
namics from the fluctuating vortex loops, to which I turn 
next. 



IV. DYNAMICS OF THE GAUGE-FIELDS 

The zero-temperature partition function for the cou- 
pled system of d-wave quasiparticles and superconduct- 
ing phase fluctuations is therefore 



Z 



D[^,S,v\e 



-(S[-f,a,v]+S u{1) [S,v]) 



(11) 



with S[^,a, v\ defined by the Eq. (10), and with 
Sjj(i) [flj v\ to be derived by integrating out the phase fluc- 
tuations. For simplicity, I will assume that these may be 
described the 2 + 1 dimensional XY model. The bare 
stiffness for the phase fluctuations will be assumed to be 
provided by the high energy modes that have been inte- 
grated out in arriving at the low-energy theory. Our goal 
will be then to rewrite the partition function for the XY 
model as the functional integral over the fields a and v. 
In particular, we want to integrate over the topological 
defects implicit in the XY model. 

I first discretize the space and the imaginary time in 
writing the partition function of the XY model. This 
is done to facilitate a more rigorous treatment of the 
topological defects, and it will prove possible to return to 
the continuum description we employed until now. On a 
lattice, in the standard lattice gauge-theory notation |l7| ] 

Z xy = / (TTd0 l )exp(A' y~] cos((f>i +(l - (pi)), 

Jo ..... 

i i,jj,—x,y 7 r 

(12) 

where the index i labels the sites of a three (2 + 1) di- 
mensional lattice, and x is the lattice unit vector in the x 
direction, for example. For simplicity, full isotropy in the 
XY model is assumed. Using the Villain approximation 
p0| and then integrating over the phases leads to 

/OO ' 1 

d5^exp(-— ^(Vxs,) 2 (13) 

i 

where = (?ij ir , n^, rii lV ) is an integer vortex-loop vec- 
tor variable, satisfying the constraint V • = (indi- 
cated with the prime on the sum). V and Vx should be 
understood as the lattice gradient and the curl, respec- 
tively. Summing over Hi forces Si to take integer values, 
and the above expression becomes the standard current 
representation of the XY model fPjv . 



Next, I imagine dividing a given configuration of vortex 
loops into two arbitrary groups, and write Hi = nA,i + 
with V ■ nA.i — V ■ ng,i = 0. We will want to 
sum over all integer Ha,% and Hb,u in order to average 
over all possible divisions of vortices into two groups. 
Introducing the lattice version of the fields Si and tfj as 
Bi + b% = 2irnA,i, Bi — bi — 2-KfiBs, where bi = V x Si 
and Bj = V x u, , I write |3l| 



d[a,v,t,s,r\ ^ exp- 

- OO -* -* 



E^( Vxs «) 2 ( 14 ) 



+i2ivsi ■ (n A ,i + H B ,i) 
+iti ■ (Bi +bi — 2-kuaa) + in(Bi — bi — 2nri Bi ))- 

The summations over nA,i and H B ,% then enforce the con- 
straints Si — U — fn A ,i, and Sj — t\ — rh B ,i, where rfiA,i 
and rriB,i are new integers. Performing the Gaussian in- 
tegrals over Si, yields 

/OO 1 
d[a,v] exp-^^Vx^i (15) 

ivi ■ (V x (rriAs + + i®i 1 (V x (tTiaa - Wb,;))]. 

This can be further simplified by noticing that the ac- 
tion is quadratic in the Volovik's field v, which can also 
be integrated out. In doing so I will neglect the addi- 
tional coupling of v to the charge current J in the Eq. 
(10), which only leads to additional irrelevant interac- 
tion between spinons. The integration over vt in the last 
equation then gives 

/OO ' ^ 

da 2J ex P"l][^ ( 16 ) 

(V x (rfiA,i + mB,i)) 2 + idi ■ (V x (rriAs - m B ,i))] 

Integrating over Si in (16) would give back the current 
representation of the XY model, Eq. (13). Alternatively, 
we can introduce the real variables and and 
write 

d[s, $_,$+] ex p-E(^( v x *+^) 2 

IaXb 

+ia, ■ (V x $_ i) + i2n(l A ,i ■ $A,i + h,i • $s,i)) 

where 3>+,-,-t = &A,i ± $b,i- The summations over 
the auxiliary link variables Ia,b force <&a and $b, and 
therefore $ + and <&_ to be integers. To preserve the 



gauge invariance ($ 



+ V M x_.i) of the last expression we must impose 
V • l A ,i = V-f B ,i = Q, Q. We may next add a small 
chemical potential for the link variables Ia,b to the ac- 
tion in Eq. (17) as the term xJ2i0A i + ^b Up to the 



G 



Villain approximation, the last expression is then equal 
to 

d[a,$ Al $ B ] / d[9 A ,6 B } (18) 



• C0S(6U,; - - 27T$A,i,c) 



L 8iT 
1 

1 

2a; 



cos(9 B ,i - 9 B ,i+o - 2-K$> B ,i,v)], 



where I introduced two sets of "dual" angles 9a,% and 6 B ,i 
to insure the gauge invariance, and imposed the "frozen" 
limit x — > 0. The integration over Si in the Eq. (18) 
together with the frozen limit ultimately sets 9a,i = $B,i, 
so the last equation becomes another representation of 
the frozen lattice superconductor (FLS), which is well 
known to be dual to the XY model in three dimensions 

In principle, one would like to integrate out all the 
fields other than a in the Eq. (18), to be left with the 
effective action Su(i) [a] for a only. The result would be 
an interacting theory for a, which can be expanded in 
powers of a, for example. Instead of doing this, I will ap- 
proximate the Su(i) [a] with the effective Gaussian action 
for a, that reproduces the gauge-field propagator in the 
full theory (18). This approximation may be understood 
as the self-consistent mean field theory for a, with the 
effect of integration over all other fields in (18) lumped 
into the form of the propagator. 

In this approximation the problem of dynamics of the 
gauge-field a reduces to the computation of the two-point 
correlation function for a from the representation of the 
XY model in Eq. (18). I therefore introduce the source 
term into the last expression by adding i ji ■ (V x 5;) 
to the exponent. Then 



((V x a) My (V x a) jtlx ) 



d 2 

— — In Z xy I j =0 . (19) 



It is convenient then to integrate over a in the Z xy first. 
One finds 

((V x a)i, v (V x a) j} u) = 8i,j8 v ,u lim (20) 

x— >0 

7T 2 

— (cos(6> 4 - 9 i+c , - 2tt§^ v )) fls , 
x 

where the last average is to be taken over the configura- 
tions of the FLS 

Z xy =\jmj d[$,%xp-^[^(V x $) 2 (21) 



(cOs(6>i - 9 l+0 - 27T$ ii £>)]. 

X 

It is well established that the lattice superconductor at 
a small but finite "temperature" x has a phase transition 



as K is varied in the same universality class as in the 
frozen limit x — (ItJ , 1 34 1 , [^6| . We may therefore relax 
the constraint x — > with impunity and assume x to 
be finite. The average that appears on the right hand 
side of the Eq. (20) can then be computed, for example, 
by using the mean-field approximation to the FLS action 
(21) (see Appendix A). This yields 

-(€05(9, - 9 l+0 - 2-K^ i , u )) FLS cx |(exp(^))| 2 . (22) 

X 

This result is quite general, and it simply expresses the 
fact that in the ordered phase of the theory (21) the dual 
angles become correlated, while at the same time the 
gauge-field becomes massive via Meissner effect. The 
gauge field fluctuations can then be neglected, which 
makes the requisite average finite when the dual angles 
9 order, i. e. in the disordered phase of the original XY 
model. 

Returning to the continuum notation, and switching 
to the Fourier space, the gauge-invariant expression for 
the correlation function (19) at low momenta is therefore 

<(V x 3)„(V x a) M ) cx (|($)| 2 + 0(q 2 ))(S^ - 

(23) 

where I allowed, in general, for some momentum depen- 
dence (the term 0(q 2 )). 0(q 2 ) term should be expected 
to appear in a more sophisticated approximation for the 
gauge-field dynamics than provided by the Eq. (20). To 
the lowest order, the integration over all other fields in 
(18) effectively yields the Maxwell term for the gauge field 
a, with the stiffness inversely proportional to the expec- 
tation value of the dual loop condensate (<&) ~ (e l6 ) that 
reflects the phase of the XY model. This is the main 
result of this section. When the dSC is phase coherent 
and the vortex loops are finite in size, ($) = 0, and a is 
infinitely stiff, and in first approximation may be consid- 
ered decoupled from spinons. When vortex loops blow 
up, (<&) ^ 0, phase coherence is lost, and the spinons are 
minimally coupled to a massless gauge field. This is in 
agreement with the physical arguments advanced in [[[si . 

At high temperature one can neglect the fluctuations 
in the imaginary time direction and deal with the purely 
2D problem of point vortices and antivortices. This sim- 
plifies the analysis in that no gauge invariance needs to be 
insured in the Eq. (17), so no dual angles are required 
p2| , Q. One then ends up with the thermodynamic 
vortex fugacity playing the role of the dual condensate 
Jl5| , and with the simpler sine-Gordon theory instead of 
the FLS. For an alternative derivation of the gauge field 
dynamics at T ^ and in continuum that is in full accord 
with the conclusions of this section I direct the reader to 
the Appendix B. 

There is an additional subtlety in going from the lat- 
tice to the continuum theory that is worth registering. 
The partition function in the Villain approximation for 
the XY model in the Eq. (17) has the symmetry under 
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o-i,fi — * Oi,/i + 27r7it )A1) with m tll integer, that becomes bro- 
ken when a small chemical potential i / for the link 
variables Ia,b (in passing to the Eq. (18)) is added. This 
periodicity would dictate that the summation over the in- 
teger vortex variables in Eq. (17) with x = should yield 
a compact term for a. Absence of the chemical potential 
x, i. e. of the vortex core energy, in the Eq. (17), on the 
other hand, must be regarded as an artifact of the Vil- 
lain representation to the original XY modeLin which, 
as well known, vortices do cost finite energy ■ This is 
because the Villain approximation reproduces correctly 
only the long-range part of the vortex interaction, while 
the short range part needs to be modified " by hand" |37]] 
in order to obtain the finite core energy. The dynamics 
of a should therefore be determined from the theory with 
x ^ 0, and by the non-compact Maxwell term, as in the 
Eq. (23). Possible effects of compactness of a on the 
picture developed in this paper are discussed in the sec. 
XII. 



V. DYNAMICAL BREAKING OF CHIRAL 
SYMMETRY 

The effective T — low-energy theory for the inter- 
acting system of d-wave quasiparticles and fluctuating 
vortex loops, after the integration over vortex loops is 
therefore 



S[9] 



d 2 fdT{*i[7o(a T + iao) + liv f {d x + ia x ) (24) 
+j2\v A \(d y + ia y )]9 1 + (l^2,x^y) 



1 



2\m 



;(c 2 (Vxot + (Vx«m 



where I omitted the higher derivative terms, and the 
terms quartic in 'P. This is the standard three dimen- 
sional quantum electrodynamics (QED3), with two im- 
portant caveats: 1) the coordinates x and y are ex- 
changed for the second Dirac field, 2) there is an in- 
herent anisotropy in the model, Vf 7^ va 7^ c, where 
c is a characteristic velocity for the phase fluctuations 
psfl . First, let us consider the simpler isotropic limit 
of the theory, Vf = va — c. There are sixteen 8x8 
matrices then that either commute or anticommute with 
the three 8x8 7-matrices that appear in the Eq. (24): 
diag{-y ,j }, diag{j!,j 2 }, diag{j 2 ,7i}- First, there are 
eight block-diagonal Hcrmitean matrices 

I®h,a3®h,I ® 735,03 ® 735, (25) 

that commute, and 

I® 13, 0-3 ®73, 1® 75,^3 75, (26) 

that anticommute with the 7-matrices. Here, 73 = a 2 ® 
o~2 , 75 = 0-3 &> I, 735 = *7375, and I4, = I ® I. Next, there 
are eight more block-off-diagonal Hermitean matrices 



1 2 
o-\ ® ^|(72 - 7i)73, 0-2 ® ^(72 - 7i)73, (27) 

°~1 ® 7^(72 - 7l)75, CT2 ® -^(72 - 7l)75, 

that commute, and 

1 1 
0-1 <8> -^=70(71 + 72), C2 <8> -^=70(71 + 72), 

that anticommute with the 7-matrices. I call these six- 
teen generators G;, i — 1, ...16, in the above order. The 
isotropic QED3 in the Eq. (24) is invariant under a 
global unitary transformation 



where 



U 



(29) 



(30) 



This follows immediately by observing that all 
the generators commute with the 8x8 matrices 

diag{j ,7o}diag{'Yi ) j2}, and diag{^ ^ }diag{-f 2 ,li\ 
by construction. The unitary transformations in (29) 
can be shown to form the Lie group U (4) . Following 
the standard terminology in the field theory literature, I 
will refer to this symmetry of the QED3 as "crural". 

As a first step towards understanding of the meaning of 
the chiral symmetry in the present context, it will prove 
useful to consider how it may be broken. QED3 is well 
known to have the chiral symmetry spontaneously broken 
ppf , by dynamical generation of the mass term in the 
action (24): 



f 2 



(31) 



withmcx |($)| 2 , i. e. proportional to the effective charge 
of the QED3. Containing just a single 7- matrix, the mass 
term in the Eq. (31) breaks all the anticommuting gen- 
erators, Gi with i = 5,6,7,8,13,14,15,16. The chiral 
symmetry is reduced from [7(4) to U{2) x 17(2), with 
eight generators preserved. The fermion mass is gener- 
ated dynamically due to the coupling to the gauge field. 
To see this, neglect the wave-function renormalization 
and the vertex corrections (which can be rationalized in 
the limit of a large number of Dirac fields N), and write 
the self-energy as 



£(g) = |<d>)| 2 7„ 



d 3 p D Ufi (p- $)E(p) 



(27T)* 



E*(p) 



7m, 



(32) 



where q — (u>,q xi q y ). The gauge-field propagator in the 
transverse (Landau) gauge is 



D Vfl (p) = (S^-pvp^/ip 2 + n(p)), 



(33) 
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where IT(p) is the self-consistently computed polariza- 
tion. At p <C S(0) = m, assuming a finite mass m gives 

D7T m 

For the polarization at all momenta see the Appendix 
C. The Eq. (32) was analyzed in [^(J (see also Appen- 
dices C and E), and there is a solution with finite m for 
the number of Dirac fields N < N c = 32/tt 2 = 3.24. Full 
numerical solution that includes the wave- function renor- 
malization and vertex corrections confirms that N c w 3 
p9[ , almost independently of the choice of vertex. Lattice 
simulations give 3 < N c < 4 ]io| , or at least that N c > 2 
KjJ . It therefore seems reasonable to conclude that for 
N = 2 the chiral symmetry in the isotropic QED$ be- 
comes spontaneously broken when the vortex loops un- 
bind and ($) ^ 0. 

Since the matrix 70 commutes with the electron-spinon 
transformation in the Eq. (9), it is easy to rewrite the 
mass term in the QED3 in terms of the original electron 
operators: 

2 

m^$A-*mr {[cU^™) c +(fc-Qi> w «) (35) 

-cL(-fc + Qi,w n )c_(-fc, uj n )] + 
[c\(k- Qi,w„)c + (fc, ui n ) 
-cL(-fc,cJ„)c_(-fc + Qi,w n )]} + (1 2). 

The reader will recognize this as the low-energy part of 
the staggered potential along the spin z-axis 

mftnr £ c„ 8 (Q,^>4(,»c, (r », (36) 

^ cr=±,i=l,2 

so the mass in the QED3 is nothing but the sponta- 
neously generated SDW order parameter. The periodic- 
ity of the SDW is set by the vectors Qi, and thus tied to 
the Fermi surface. The SDW order is established as soon 
as the phase coherence is lost, and the charge ($) 7^ 0. 
In the large-N approximation [|2(J one finds that 

m m 16| (*) fe-^V^W^V. (37) 

Since N c w 3, for A = 2 one finds that m - 1(T 2 |($)| 2 . 
This extreme "lightness" of fermions in the QED3 de- 
rives from the fact that the mass comes solely from the 
interaction with the soft gauge-field. 

Breaking of chiral symmetry in the QED3 also implies 
that the energies of spinons have become complex and 
finite in the phase incoherent state with ($) 7^ 0. In the 
simplest approximation the electron propagator may be 
computed as a product of the spinon and the gauge-field 
propagators, so a spinon "gap" should imply a charge 
gap as well, i. e. the system becomes an insulator [ p2[ . 
In section IX I discuss how spinons should actually be 



confined in the insulating state. Staggered magnetiza- 
tion, charge gap, and the spinon confinement when taken 
together imply that the state with broken chiral symme- 
try is nothing but the standard, albeit a weak, SDW. It 
seems natural to assume then that this state is contin- 
uously connected to the antiferromagnet near half-filling 
in cuprates. This expectation is further corroborated by 
considering the effect of Coulomb interactions, which is 
done in section VIII. 

It has been already mentioned that we have some free- 
dom in choosing the representation of the 7-matrices. In 
particular, it was the specific choice of 70 that led to the 
cos-SDW order parameter displayed in the Eq. (35-36). 
In the next section I discuss how "rotating" the cos-SDW 
by the broken chiral generators leads to a different insu- 
lating states. 



VI. MORE ON CHIRAL SYMMETRY: THE 
SPACE OF INSULATORS 

In discussing the pattern of chiral symmetry breaking 
in the QED3 one needs to distinguish at least two dif- 
ferent cases. The isotropic theory (va = v f) has the full 
U(A) symmetry in its massless phase, so the mass term 
breaks eight of its sixteen generators. In cuprates fl43| , 
however, Vt/v^ ~ 10, and even with m = the symme- 
try is only U(2) x U(2), generated by the block-diagonal 
Gi i — 1,...8. How such a large anisotropy affects the 
value of N c is a non-trivial problem, and is addressed in 
a separate publication jwj] . Here I will consider only the 
effect of anisotropy on the chiral symmetry, and assume 
it is reduced to U(2) x U(2). It suffices then to look at 
each Dirac component in the QED3 separately, i. e. con- 
sider just the 4x4 representation of the 7-matrices, as 
defined right below the Eq. (5). 

It can be easily shown that any matrix that anticom- 
mutes with both M\ and M% and squares to unit matrix 
may be chosen as 70, and will lead to a representation 
of the 7-matrices like in the Eq. (6). The mass term 
~ m^*^/^ in the action would then gap the quasiparti- 
cles, in analogy to the standard relativistic Dirac equa- 
tion. The problem of different chiral orders is therefore 
nothing else but finding all the ways in which d-wave 
quasiparticles can spontaneously acquire such a "rela- 
tivistic mass" . It will be useful to introduce " directions" 
in the space of broken symmetry states, as a set of lin- 
early independent matrices that anticommute with M\ 
and M2, and square to one. It is easy to show that in the 
4x4 representation there are only four such matrices 

7o,73,75,«7i72, ( 38 ) 

with 70 = 70, and where 71 = — iM\, 72 = 1M2, 73 = 
03 ® 02, and 75 = (T2 <8> /. In principle, any of these four 
if used instead of our 70 in the construction of the Dirac 
theory in the Eq. (6) and in the mass term would give a 
relativistic gap to Dirac fermions. The last matrix, 
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(39) 



however, being a product of two 7-matrices does not 
break the chiral symmetry, and is believed not to be 
spontaneously generated in the QED3 p0| ], |45| ], I there- 
fore focus on the remaining three. Choosing one among 
{70,73,75} as the 70-matrix in the mass term reduces 
the SU C (2) subgroup of U{2){= U(l) x SU C (2)), gener- 
ated by {73,75,735}, to U c {\). The two anticommuting 
generators of the SU C (2) that are broken then rotate the 
chosen order parameter towards the two remaining "di- 
rections" in the chiral space. For example: for our choice 
of 70 =7o, it is 735 that remains unbroken in the cos- 
SDW phase, whereas the broken generators rotate the 
cos-SDW order parameter as 

e «07^ oe -i0 7i = cos(26»)7o - sin(26')7 i ;i = 3, 5. (40) 

Choosing i — 5, for example, for both Dirac fields rotates 
the cos-SDW in the Eq. (35) into 

m I d 2 fdr cos(Qi ■ f + 2Q)ac\(f, r)c cr (r, r). 

J a=±,i=l,2 

(41) 

Chiral rotations generated by 75 thus correspond to slid- 
ing modes of the SDW. 73, on the other hand, ro- 
tates 70 towards the direction of 73, which describes an 
additional particle-particle pairing between the neutral 
spinons, with the opposite sign for the diagonally op- 
posed nodes. This may be understood as an additional p- 
wave pairing between the spinons, so the state described 
by 73 order parameter may be called the "d+ip" state 
p6| . This state preserves the superconducting U(l) sym- 
metry and the translational invariance, but breaks the 
spin-rotational invariance and is odd under parity. Since 
73 does not commute with the electron-spinon transfor- 
mation (9), however, "d+ip" state can not be that simply 
expressed in terms of electronic operators, as it proved 
possible for the SDW states. The relationship between 
the directions in the order parameter space {70,73,75}, 
and the chiral generators may be summarized pictorially 
as on Fig. 3. 

It is instructive to look more closely at the origin of 
the U(2) symmetry (per Dirac component) that appears 
in the low energy theory of the dSC. First, the transfor- 
mations in the U{1) subgroup of U{2) = U(l) x SU C {2) 
are analogous to the spin rotations around the z-axis. To 
see this, consider the conserved current that corresponds 
to the U(l) subgroup: Ji.^ — ^ij^i, so that the con- 
served charge is simply the z-component of the total spin 
of the low-energy quasiparticles, one charge per each pair 
of nodes. Of course, the quasiparticle action (2) has the 
full 5*0(3) spin symmetry, and this is not to imply that 
a part of it is broken in the dSC. It only means that in 
writing the full action (2) in terms of the Dirac fields (3) 
only the subgroup of spin rotations around z-axis is rep- 
resented by a simple 4x4 matrices that act on "J. The 




FIG. 3. The corners of the triangle represent the three chi- 
ral directions in the space of insulating states that descend 
from the d-wave superconductor. At the side opposite to 
a particular direction lies the corresponding unbroken chiral 
generator, while the remaining two broken generators rotate 
the chosen insulator towards two other directions. 

rest is still present, but not that obvious in our choice 
of the Dirac fields, which was made to make the chiral 
symmetry manifest. (For a complementary representa- 
tion that is fully rotationally symmetric at the expense 
of chiral symmetry see the Appendix D.) The U(l) sub- 
symmetry is therefore always present, both in the super- 
conducting and the insulating states. The SU C (2) factor 
is more interesting. The conserved currents (per pair of 
nodes) in the dSC that correspond to this symmetry are 
■A> = ^il^^i, where T = 73,75,735. As we have seen 
already, the 75 generator simply translates in the diago- 
nal direction. The corresponding conserved charge may 
be written as 

Qi = J drdTjf Q =T ^ ±c\(k,ui n )c <J (k,0J n ), 

(42) 

and may be identified with the quasiparticle momentum 
along Ki. More precisely, under the translation of the 
original electron operators c CT (fc,(jj) — > e lk ' R c a (k,oj), the 
spinon field transforms as 

* 4 (r,r) -> e^K^iir + R,t), (43) 

where k = Ki + q. The low-energy theory therefore has 
more symmetry than the original action (2), as the chi- 
ral rotation by 75 and the translations of \& separately 
are still the symmetries of the theory (24), while only 
when combined as above do they represent an ordinary 
translation. Nevertheless, breaking of chiral generator 75 
implies breaking of the translational symmetry in the the- 
ory. The remaining two generators of the chiral SU C (2), 
73 and 735, on the other hand, are not related to any spa- 
tial symmetry. They should be understood as " internal" , 
and approximate, symmetries of the dSC that emerge at 
low energies. They rotate the translationally invariant 
"d+ip" state into a SDW, and therefore connect the two 
fundamentally different types of insulators. 



10 



The reader should also note that in the action (24) the 
order parameter can be rotated independently for the 
first and the second Dirac field. Any linear combination 
of 70 73, and 75 is a regular order parameter too. Since 
75 is just a sin-SDW, the fundamentally different states 
are just the SDW are d+ip state. This, however, already 
leads to a variety of insulating phases. For example, one 
can choose the cos-SDW for the first Dirac field (Qi), 
while being in the d + ip state for the second. This would 
correspond to a one-dimensional SDW along one of the 
diagonals. 

With the velocity anisotropy neglected the QED 3 has 
a larger U(4) symmetry, with sixteen generators G{. The 
mass term now breaks all eight anticommuting gener- 
ators, and the chiral manifold of insulating states be- 
comes larger. For instance, rotating the cos-SDW with 
9 = 7r/4 and the generator G = G13 — G15 leads to a 
uniform state with an additional " s" component of pair- 
ing between spinons, "d+is" [Q. Interestingly, rotating 
the cos-SDW by block-off-diagonal generators may also 
lead to charge stripes. For example, taking 9 = 7r/4 and 
G15, rotates the 8x8 cos-SDW order parameter / <g> 70 
into (l/2)<7i (g) (71 + 72). When written in terms of the 
electronic operators, this order parameter corresponds to 
the one- dimensional charge density- wave with the peri- 
odicity Pb = K2 + K\, and with residual pairing corre- 
lations in the orthogonal a-axis direction. It has been 
known that stripes indeed occur in some high-T c mate- 
rials [Q . Here they emerge as insulating cousins of the 
d-wave state in the isotropic limit of the theory. It is also 
interesting that stripes seem always to be accompanied 
by the residual pairing correlations, so one can think of 
them as weakly coupled one-dimensional systems on the 
verge of becoming phase coherent. 



VII. REDUCTION OF CHIRAL SYMMETRY BY 
THE IRRELEVANT TERMS 

We saw that the low-energy theory of dSC has the chi- 
ral U{2) symmetry per Dirac component, which when 
spontaneously broken leads to emergence of the SDW or 
the d + ip insulators. This enlarged symmetry arises only 
at low-energies, and the irrelevant terms omitted in the 
Eq. 6 reduce the U(2) to U(l) x U c (l). In this section I 
show the higher order derivatives and the Hubbard repul- 
sion reduce the chiral SU C (2) symmetry to just transla- 
tions, generated by 75. However, we will also find that if 
both perturbations are weak it will actually be the SDW 
solution that is energetically preferred. 

Let us first consider the higher derivative terms in the 
Eq^ (6). Since £{k - Q x ) = £{q - K x ), and £(q - K x ) = 
£(i?i — q), and analogously for A(fc), one can write the 
second-order derivatives in the Eq. (6) as 

Si = -i f d 2 rdT*' l75 ( 7i r(a 2 ) - 72 A"(d 2 ))*i (44) 



+(1 -» 2,x~y), 

where £" and A" are the functions coming from the 
expansion of £(fc) and A(fc) around Ku2), respectively. 
Their specific forms are model dependent, and will not be 
of importance here. What is important is that S\ man- 
ifestly breaks the part of the chiral symmetry generated 
by 73 and 735 , while preserving only the translational in- 
variance, generated by 75. One can easily prove that the 
same holds to all orders in the gradient expansion. 

Next, consider the Hubbard-like short-range repulsion 
term, in the continuum notation, 

H V = U J d 2 x{n+{x) + n_(f)) 2 , (45) 

with U > 0. Retaining again only the excitations near 
the four nodes, one can write this as 

Su = U I d 2 rdr[i ^757i^] 2 - (46) 

■* 4=1,2 

The reader is probably not surprised that it is again only 
75 that remains the symmetry generator. This is because 
75 in our formalism is related to translations, which are 
always the exact symmetry of the action (2). So both the 
higher derivative terms, and the quartic repulsion term 
reduce the chiral SU C (2) subgroup of U(2) to U c (l), the 
translations. One could therefore naively expect that it is 
the d+ip state, which is translationally symmetric, that 
may be preferred by these perturbations. To decide on 
this, however, it is not enough just to know the symmetry 
of the action, since the new terms Si and Su may turn 
out to disfavor the d+ip state. Assuming that both £" 
and A" terms are small, one finds that the contribution 
to the energy of the SDW (or d+ip) state is of second 
order in the S\. The interaction term, on the other hand, 
yields 

i 

with the average taken over the massive QED3 with 
7o = v\ ® I (cos-SDW). The result, of course, is the 
same for the sin-SDW, or for any linear combination of 
the cos-SDW and the sin-SDW. Alternatively, if one as- 
sumes the d+ip ordering, one finds that Su gives then a 
positive contribution to its energy, to the first order in U. 
Although Sjj is only translationally symmetric, it actu- 
ally inhibits the formation of the translationally invariant 
state, and prefers the ordering to be in the " orthogonal" 
direction, i. e. the SDW. 

If both the interaction and the gradient terms are 
weak, it will therefore always be the SDW solution that 
is energetically preferred. This is because both the repul- 
sive and the higher derivative terms are equally irrelevant 
by power counting (and have the engineering dimension 
— 1), the gain in energy due to SDW is of first order only 
in U. The gradient terms affect the energy of the SDW 
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only to the next order. So at long enough length scales 
one can alway neglect higher-derivative terms as com- 
pared to the repulsive interaction, which then serves to 
select the SDW over the d+ip insulator. 



VIII. MEAN-FIELD THEORY WITH REPULSIVE 
INTERACTION 

The message from the previous section is that the quar- 
tic term that represents a short-range repulsion, although 
irrelevant, at low but finite energies is still finite, and it 
breaks the chiral symmetry in favor of the SDW state. 
This is its first important role. The second is that once 
the chiral symmetry is dynamically broken by unbind- 
ing of vortex loops, the quartic term affects the size of 
the order parameter, and therefore sets the scale for the 
value of the SDW transition temperature. In this section 
I formulate the simplest mean-field theory of the chiral 
symmetry breaking in presence of the repulsion term, and 
demonstrate that it drastically increases the value of the 
SDW order parameter at T = 0. 

We have seen that unbinding of vortex loops leads 
to weak SDW order, but with the order parameter or- 
ders of magnitude smaller than the coupling constant 
|($)| 2 . Assuming that the dual condensate as a function 
of doping x should be of the same order of magnitude 
as the superfluid density on the other side of the tran- 
sition (at x — x u ), \(3>(x u — S))\ 2 ~ p s f(x u + <5), and 
that Uemura scaling Q T c (x) cx p s f(x) is obeyed, the 
identification of the size of the SDW order parameter 
with the transition temperature T sc i w (x) suggests that 
T s dw{x u — 8) <ti T c (x u + 5). The difference in the rel- 
evant scales for the superconducting and the SDW or- 
derings is in accord with the known phase diagram in 
the underdoped regime. Starting from half-filling, with 
increased doping the antiferromagnetic order is quickly 
lost, and only at a larger doping the dSC appears. I at- 
tribute the absence of the obvious SDW order very near 
the superconducting phase to the inherent weakness of 
the spontaneous chiral symmetry breaking in dSC. As- 
suming that the weak SDW smoothly evolves into the 
commensurate antiferromagnet near half-filling, the ob- 
vious problem then becomes the following: how should 
one understand the dramatic increase of T sc i w (x) near 
half filling, all the way up to ~ 300K? 

The answer is provided by the observation that al- 
though the repulsion U is irrelevant if weak enough, it 
enhances the SDW order once it became spontaneously 
generated trough the interaction with the gauge field. 
To show this I will consider the mean-field theory of the 
QED% with the additional Su quartic term. First, notice 
that in the Hartree-Fock approximation the Su term gets 
replaced by the effective quadratic term 



Su 



U{%%) / d 2 rarTr(f^7 57l 757i^): (48) 



FIG. 4. The SDW order parameter m in units of 
|($}| 2 as a function of dimensionless short-range repulsion 
5 = [/|<$>| 2 /(2^) 2 . 



with the average to be calculated self-consistently within 
the theory quadratic in fermionic fields. The above term 
corresponds to the decoupling in the exchange (Fock) 
channel, since the direct (Hartree) term vanishes. There- 
fore in the Hartree-Fock approximation, after the Franz- 
Tesanovic transformation 



(49) 



Assuming a uniform x = ~U \) o , and treating 
the gauge-field fluctuations in the large-N approximation 
leads to two coupled equations for x an d for the momen- 
tum dependent fermion self-energy 



x = u 



d 3 k 



E(fc) 



(2tt) 3 k 2 + E 2 (fc) : 



= x + J 



d 3 k 



2|($)| 2 E(fc) 



(2tt) 3 (fe 2 + £ 2 (fc))(p 2 + n(p)) 



(50) 



(51) 



with p = k — q. When [7 = these reduce to the Eq. (32), 
which leads to N c = 32 /ir 2 . When loops are bound and 
($) = 0, on the other hand, E(g) = x, and the Eq. (50) 
allows a nontrivial solution only when the dimensionless 
coupling g = UA/(2n) 2 > 1, where A is the UV cutoff, 
A < T*. Assuming that long-range SDW order and dSC 
do not coexist, I take that g < 1 in the superconducting 
phase, so that the quartic coupling is there irrelevant. 
With (4>) ^ 0, however, small g ceases to be irrelevant, 
since there is now a small mass scale to effectively cut off 
its flow. Since £(<?) is quickly damped for q >> K^)! 2 , 
one can take the UV cutoff in the above equations to be 
A~ K^)! 2 - The above equations were studied before @ , 
|pl[ in the context of gauged Nambu-Jona Lasinio model 
of chiral symmetry breaking in particle physics. Here I 
solve the equations numerically for N = 2, as discussed 
in the Appendix E. The result is presented in Fig. 4. The 
main point is that as the superconducting phase is more 
and more disordered and the dual condensate grows, the 
presence of a moderate repulsion between electrons in- 
creases the SDW order parameter at T = by one to 
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two orders of magnitude. Recalling the above argument 
that compares T$dw to the superconducting T c on the 
other side of the superconductor- insulator transition, this 
appears to be in qualitative agreement with the generic 
behavior observed in underdoped cuprates. 



IX. CONFINEMENT OF SPINONS 

In the superconducting state, the electrically neutral 
low-energy spinons represented by the fermionic field 'J 
in the QED3 are well-defined excitations. This effective 
spin-charge separation implicit in the superconducting 
state was emphasized in j52[ , and more recently in [Q. 
One may therefore naturally wonder if this form of spin- 
charge separation will survive once the superconductiv- 
ity is lost via unbinding of vortex loops. The answer 
seems to be no. It is believed that the chiral symm etry 
breaking and confinement go together in the QED 3 [g0| . 
Qualitative argument why it should be so is provided by 
the low momentum form of the polarization tensor in 
the Eq. (33) ||, §4) : U(q) ~ q 2 /m for q < m, so 
in two dimensions spinons are at large distances bound 
by a logarithmic potential. One may independently ar- 
rive at the same conclusion by analytically continuing the 
fermion propagator in the broken symmetry phase to real 
frequencies |5^] , to find that its poles lie at complex ener- 
gies with both real and imaginary parts proportional to 
the chiral mass. The chiral symmetry breaking and con- 
finement of spinons seem therefore to go hand in hand 
in the QED3, so the states with broken chiral symmetry, 
including most importantly the SDW, should not have 
well defined fermionic excitations even above the mass 
"gap". 

Dissapearance of spinons from the spectrum in the 
insulating phase, if required, imposes a rather non- 
trivial constraint on a candidate theory for underdoped 
cuprates. For example, one could imagine a completely 
different mechanism of chiral symmetry breaking in dSC: 
even without the gauge field, simply increasing the quar- 
tic coupling U above a certain value (U c A/(2tt) 2 = 1 in 
the Hartree-Fock approximation) would open the gap for 
spinons and lead to SDW order. This would be analo- 
gous to the chiral symmetry breaking in the Nambu- Jona 
Lasinio and related models @, @, @, H- The cru- 
cial difference, however, is that such a mechanism would 
yield well defined spinon excitations at energies above 
the gap, in the insulating state. The integrity of the 
gaped spinons is assured essentially by the Landau phase 
space arguments. Such a "deconfined" antiferromagnet 



was dubbed AF* and studied in 24 1, for example. From 
this point of view it becomes a non-trivial problem to un- 
derstand how spinons could be removed from the spec- 
trum. In the QED3 this is accomplished via the same 
non-perturbative mechanism that yields chiral symmetry 
breaking, described by the Eq. (32), for example. 

Having said all this, it needs to be realized that in 



a weak SDW confinement of spinons is effective only 
over very large distances, L » 1/m. At intermedi- 
ate scales, the polarization H(q) ~ q, so the potential 
between spinons is ~ 1/r, and at intermediate distances 
1/m >> L >> 1/1 ($)| 2 spinons will appear effectively 
deconfined. In this sense it is still meaningful to think 
about underdoped cuprates as exhibiting an effective 
spin-charge separation. Computing the electron spectral 
function by taking the gauge-field fluctuations into ac- 
count in large-N approximation fl§| l, which suppresses 
the dynamical symmetry breaking, for example, gives re- 
sults in qualitative agreement with the experiment J60| . 
As one continues to underdope, however, SDW order pa- 
rameter grows, and spinons become more strongly con- 
fined. In the strong antiferromagnet at half-filling there- 
fore, on may expect spinons to be confined already at 
atomic distances. 



X. EXPERIMENT 

The principal consequence of the QED3 theory of un- 
derdoped cuprates is, of course, the antiferromagnetism 
itself. All the materials that become d-wave supercon- 
ductors with doping are insulating antiferromagnets in 
its parent state. Furthermore, the sharp spectral fea- 
tures in the dSC should become very broad in the in- 
sulator, since there is a soft (propagator ~ l/<? 2 ) gauge 
field in the problem. Nevertheless, an insulator that de- 
rives from a dSC should partially inherit the d-wave form 
for its "gap", except for its finite value in the nodal 
directions. This is in very good agreement with the 
ARPES measurement on the insulating Ca2CuOiCli^ 
and Sr 2 Cu02Cl2, in its parent state ||, |J. In Fig. 
5 I compare the ARPES data for the "gap" measured 
from the top of the lower Hubbard band in the insu- 
lating state with the simplest functional form consis- 
tent with the chiral mass: at the remnant Fermi surface 
w = ((E max (cos(k x ) - cos(fc„))/2) 2 + E^J 1 / 2 , where 
the chiral mass m = E mm = 75meV is chosen to be the 
T = sublattice magnetization for J = 125meV. Best 
fit is obtained then for E max = 420meV. The quality of 
the fit is actually not very sensitive to some variations in 
E min and the corresponding E max . 

The key prediction of this work is that the above 
" gaped d-wave" form of the insulating " gap" is a generic 
feature of the insulating state. Upon underdoping 
the ARPES should show the standard d-wave gap for 
sharp quasiparticlcs in the superconducting phase, which 
should evolve into a gaped d-wave form for broad ARPES 
shape in the insulating state, with the "gap" increasing 
as one approaches half filling. The rounding of the data 
at low energies should therefore be intrinsic to the insu- 
lating state, and should weaken with doping. Although 
the initial experiment on CaiCuOiCli only indicated 
such rounding, later measurements on SY2CUO2C/2 with 
higher resolution pi clearly showed the deviation from 
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FIG. 5. ARPES results for Ca 2 Cu0 2 Cl 2 (bars) and 
SriCuOiCh (dots) with E = E(k) - E(n/2,n/2) in meV. 
The line is the function described in the text. 



the simple d-wave cusp at lowest energy. More recent 
measurements fjlf indicate that the rounding of the data 
at low energies is a robust feature. It would clearly be 
desirable to perform a systematic study of this effect at 
variable doping. 

It may also be worth mentioning that some signs of 
the gap rounding in the insulator may be observable al- 
ready in the superconducting state. In BiTlYl |32|, for 
example, as one underdopes, the d-wave gap continues to 
show the cusp at zero energy, but with the slope (velocity 
va) decreasing, in spite of the increase in the overall gap 
magnitude in (tt, tt) direction. It is tempting to interpret 
this effect as a precursor of the dynamical mass gener- 
ation. Detailed study of this effect and of the spectral 
features in the insulator is deferred to a future work. 



XI. CONCLUSION AND DISCUSSION 

In summary, I have shown that the minimal theory 
that describes unbinding of vortex defects in the d-wave 
superconductor at T — is the two-component, 2+1 di- 
mensional QED, with the vortex condensate playing the 
role of " charge" . With the loss of phase coherence, the 
d-wave superconductor suffers the spontaneous breaking 
of the low-energy "chiral" symmetry, which results in a 
weak SDW order. It was argued that with underdop- 
ing this SDW smoothly evolves into the strong antifcr- 
romagnet near half-filling, with the selection and the in- 
crease of the SDW order parameter being provided by 
the repulsion between electrons. I argued that spinons 
are marginally confined in a weak SDW, and may appear 
effectively deconfined over intermediate length scales in 
the pseudogap regime. Finally, it was proposed that 
the rounded d-wave form of the "gap" in the insulating 
CazCuOiCli observed by Ronning et al. may be a con- 
sequence of the chiral mass for the approximate spinon 
excitations, as implied by the QED3. 

The present theory is similar in spirit to the approaches 
of refs. [Jt| , and |24|] , in that it attempts to understand the 
phase diagram of underdoped high temperature super- 



conductors beginning from the superconducting phase. 
It differs, however, in its conclusions to what the ground 
state that results from unbinding of topological defects 
in the d-wave state is. Whereas it was argued in pj and 
p3 that the relevant description of this process is pro- 
vided by the Ising (Z2) gauge theory, and that the re- 
sulting state may show spin-charge separation, I argued 
that unbinding of defects of unit vorticity leads to the 
dynamical symmetry breaking in the low-energy theory, 
and the accompanying confinement of spinons in the in- 
sulating state. In fact, if one demands that the insulat- 
ing state near half-filling is the standard antifcrromagnct 
with spin-1 excitations and confined spinons, the form of 
a single theory that would be able to describe both the 
dSC and the insulator becomes severely restricted. The 
QED 3 in this paper is one such theory. 

A variation of the QED3 as an effective theory for un- 
derdoped cuprates has also been considered before p3[ , 
p4j| , [ |66| as the theory of low-energy fluctuations 

around the 7r-flux phase in the large-N version of the 
Heisenberg model. In that approach the gauge invari- 
ance reflects the local particle number conservation at 
half-filling, and the gauge-field has no dynamics on its 
own. As a result, the gauge field is necessarily com- 
pact, and the theory is infinitely strongly coupled. Not 
much is definitely known about such a lattice gauge the- 
ory, which greatly diminishes its utility. Nevertheless, 
it was argued that neglecting the instanton configura- 
tions would restore the antiferromagnetic order at half 
filling, via spontaneous breaking of a different "chiral" 
symmetry, which in this case is actually an enlarged spin 
rotational symmetry [[34j, [p5|j . While this logic may at 
first appear close to the one in the present work, there 
are crucial differences. First, I begin from the supercon- 
ducting state, away from half-filling, with the gauge field 
describing vortex fluctuations. As a result, the gauge 
field is weakly coupled and non-compact near the dSC- 
SDW transition. Also, the SDW phase that obtains from 
chiral symmetry breaking may be incommensurate, and 
the approximate chiral symmetry of the low-energy the- 
ory is unrelated to spin rotations. 

Nevertheless, it may be possible to understand the 
QED3 as a low-energy description of the microscopic t-J 
model of cuprates. Starting from the mean-field slave- 
boson theory of the t-J model and integrating the con- 
straints of no double occupancy, for example, leads to an 
effective theory of the form quite similar to the QED3 
|37j, but with the Volovik's field v only. Including vor- 
tices would then be expected to introduce the gauge field 
a, as shown in this paper. The point is that irrespec- 
tively from the underlying microscopic model the theory 
of the fluctuating dSC should assume the QED3 form. 
Values of the parameters, however, may strongly depend 
on the microscopic physics: the bare stiffness K in the 
XY model for the phase fluctuations (Eq. (12)), for ex- 
ample, should be proportional to doping x in the doped 
Mott insulator [Q. Also, the charge of quasiparticles 
(the coefficient in the last term in the Eq. (10)) would 
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be expected to change from unity to ~ x, at small dop- 
ings. 

There exist further parallels between the QED3 and 
the gauge theory of the t-J model. One may formulate a 
representation of the t-J model with a U(l) gauge field 
that minimally couples to spinons and holons. It was ar- 
gued |6j| that the effect of holons would be to screen the 
temporal component of the gauge field, which then may 
be shown to halve the critical number of spinon species 
for the chiral instability, N c — > N c /2. That way one 
could avoid the chiral transition at N = 2 (assuming that 
N c w 3), and have a spin- liquid as the ground state in the 
underdoped regime instead. The tacit assumption, how- 
ever, is that uncondensed bosons (holons) at T = may 
exist in a compressible state. If the system becomes insu- 
lating with the loss of phase coherence, however, bosons 
would become incompressible and the above argument 
breaks down. This is indeed the case in the QED3: with 
the proliferation of vortices the system becomes insulat- 
ing, and all the components of the gauge field become 
massles. The same conclusion would be reached within 
the gauge theory of the t-J model if one would consider 
the incompressible state of slave bosons |38| . 

The present work shares the same philosophy with the 
recent one Q, p6| , where the massless U(l) gauge field 
as an effective description of unbound vortex loops was 
also considered. While the authors jlq] considered the 
large-N limit of the QED3, and thus precluded the chiral 
symmetry breaking, my main point is that at T = the 
spontaneous formation of the chiral condensate is noth- 
ing else but the SDW instability of the d-wave super- 
conductor. The results of the ref. |Q may therefore be 
understood as applying to the finite-T phase much below 
the pseudogap scale T* in Fig. 1. 

The problem of phase disordering of dSC has also been 
recently studied by Ye |6{|. Working in the Anderson 
gauge |28|] in which 4> s a = 4> s , 4>sB = in the Eq. (8), the 
author concluded that the gauge-field a is always massive 
when charge fluctuations are included. It is easy to see 
that this is a direct consequence of the gauge choice: in 
the Anderson gauge a — v, and not only v, but a too is 
ultimately coupled to the charge current. In my gauge 
invariant approach, on the other hand, a is completely 
decoupled from charge, and couples only to spin. In- 
clusion of charge fluctuations therefore does not make a 
massive, but simply adds an irrelevant quartic coupling 
to the QED3 Lagrangian. 

The intimate relationship between the d-wave super- 
conductivity and antiferromagnetism is also the main 
theme of the SO(5) theory of Zhang |70t| . The present 
work echoes some of that general idea, but is based on 
entirely different physical principles. In particular, al- 
though there should be a direct dSC-SDW transition in 
the phase diagram, this appears unrelated to the SO (5) 
symmetry, but comes as a consequence of the chiral sym- 
metry that emerges at low energies in the d-wave su- 
perconducting state. It is the spontaneous breaking of 
this hidden approximate symmetry that implies then the 



breaking of the spin rotational symmetry in the SDW 
phase. 

Marginal confinement of spinons we found in the weak 
SDW phase is very much in line with the speculations 
of Laughlin JTl ] , |7^] on parallels between the antiferro- 
magnetism and the confinement in strong interactions. 
In fact, the QED3 shows precisely how chiral symmetry 
breaking, i. e. SDW ordering, binds spinons into spin- 
1 objects. Deconfinement in this theory seems indeed 
tantamount to the absence of chiral symmetry breaking. 
In this context, it may be interesting to note that the 
d + id state, that would correspond to the 17172 matrix 
in (38), could lead to deconfined spinons. This state is 
outside of the chiral manifold, and it is believed that it is 
not spontaneously induced in the QED3 [^o| , because of 
the Chern-Simons term that becomes generated for the 
gauge-field. With the Chern-Simons term, on the other 
hand, the gauge-field propagator behaves like ~ q at low 
momenta, and thus spinons may become deconfined Jr3) . 
The chiral symmetry breaking in the QED3 is therefore 
nothing by the effective description of the spinon confine- 
ment. 

It is also interesting to note that were the critical num- 
ber of fermions N c < 2, the result of phase disordering 
of dSC would be quite different. Instead of symmetry 
breaking and confinement one would find a gapless, chi- 
rally symmetric state, in which spinons would be decon- 
fined. This is again because the polarization tensor would 
then be ~ q at low momenta, i. e. the interaction be- 
tween spinons would be ~ 1/r at large distances. This 
state would be similar in spirit to the "nodal liquid" W \, 
or analogous to the "algebraic Fermi liquids" fig] , |p0 |, 
]74j proposed in literature as candidates for the pseudo- 
gap phase. It has been proposed recently that N c = 3/2 
exactly ]75| , although all the actual calculations based 
on Schwinger-Dyson formalism lead to N c > 3. If N c 
is indeed that low, phase disordering of the dSC would 
first lead to the deconfined pseudogap phase, which only 
later would turn into the confined SDW phase, presum- 
ably due to the repulsive quartic term which is know to 
increase N c []50| , |pl[ . At this time it is hard to say which 
one of these two scenarios is realized in cuprates. 

The main point made in this paper is that unbinding 
of vortex loops in a d-wave superconductor at T = 
results in SDW order. It then appears natural to assume 
that the cores of fluctuating vortices are already in the 
insulating state. This speculation is in accord with the 
recent STM, neutron scattering, and NMR experiments 
H, 0> & JM, the SO(5) proposal @, |o), and 
the mean-field |81f| and the finite size QED3 calculations 
flS2| . The superconductor-insulator transition would then 
be the result of the decrease of the bare stiffness K in the 
XY model with underdoping, since K ~ x in the doped 
Mott insulator [|S7|. 
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XII. FURTHER PROBLEMS 



XIII. ACKNOWLEDGEMENT 



I finish with a tentative list of problems opened by this 
work. 

1) The role of strong anisotropy v//va 3> 1 that exists 
in cuprates is unclear. In particular, since anisotropy 
on the bare level is marginal, it may affect the value 
of N c . The preliminary results, indicate, however, that 
weak anisotropy is irrelevant, so that one would expect 
N c to be unaffected by it Q . 

2) Nature of the various phase transitions in the the- 
ory is also of interest. Whereas one expects that gap- 
less quasiparticles do not change the Kosterlitz-Thouless 
universality class of the finite temperature superconduct- 
ing transition, the nature of chiral symmetry breaking 
at finite temperature and its possible interplay with Neel 
transition is far less clear |S3| ]. In particular, in relation to 
the Uemura scaling [[l9| , one would like to understand the 
quantum superconductor-insulator criticality and how it 
may be affected by gapless spinons. 

3) Can long-range SDW and SC order coexist? In the 
approximation employed in the present work, the gauge- 
field a is considered decoupled from spinons in the dSC 
phase. This is likely to underestimate the effect of a, and 
a better approximation for the gauge-field propagator is 
needed to study its effect inside the dSC. This could be 
important in light of the recent experimental data J78[ , 
[ B4J that may be interpreted as indicating the coexistence 
of the SDW and the SC orders in some compounds p5[ |. 

4) The present work also points to a new route towards 
a deconfined phase in two dimensions: lowering N c be- 
low two would allow for an insulating phase with decon- 
fined spinons. At present, however, it is not clear how 
to achieve this within the QED3, unless the Schwinger- 
Dyson equations systematically overestimate N c |75| . 

5) The computation of the electron propagator within 
the QED3 is an important problem |42|| . This would be 
necessary for a detailed comparison of the theory with 
the ARPES measurements. 

6) As mentioned at the end of sec. IV, in the Villain 
approximation to the XY model, the gauge-field a ap- 
pears to be compact, in contrast to the Volovik's field v. 
Although this should be an artifact of the Villain approx- 
imation, it would still be interesting to understand the 
effect of compactness of a on the chiral symmetry break- 
ing in the QED3. It has been argued that the coupling to 
gapless spinons makes the single instanton anti-instanton 
pair that derives from compactness of a bound above the 
certain number of spinon components Ni nst H], p6| . 
Ni ns t may be made smaller than N c for chiral symme- 
try breaking by a large anisotropy j74j, for example. It 
is unclear, however, whether this conclusion survives the 
effects of screening by other pairs Wh - Also, even if the 
instantons can be made irrelevant above N c , below N c 
one would expect them to become relevant again with 
the opening of the spinon " gap" . This in turn could have 
profound consequences for the spinon confinement. 
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XIV. APPENDIX A 

I present the self-consistent mean-field theory of the 
lattice superconductor (20) J34|, and use it to approxi- 
mately compute the correlator appearing in the Eq. (20). 
By the Bogoliubov inequality: 

Z xv >Z Q e-^ H - H ^\ (52) 

where Z xy is the partition function in the dual form (21) 
with a finite "inverse temperature" x, and the average in 
the exponent is performed over a local mean-field Hamil- 
tonian 



4Kn 2 ^ 



\ " <j> 2 . 



(53) 



The optimal values of the parameters h and m that maxi- 
mize the right hand side in the Bogoliubov inequality are 
then determined by the equations: 



6 A h(h) 
x I (h)' 



2 K^h(h) 
m = , 
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3 h{h) 



d 3 k 



exp[ 
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(2tt) 3 F{k) + m 2 



(54) 



(55) 
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where F(k) — J2v( elk " — l) 2 ' an d t ne integral over k is 
taken over (— 7r, it). Iq and I± are the Bessel functions. 
These equations can be solved graphically, and describe a 
discontinuous transition from the phase with h = m = 
(bound vortex loops), to the condensed phase ^ 0, m ^ 
(infinitely large vortex loops) . 

The requisite average in the Eq. (20) is easy to com- 
pute in the mean-field theory that has different sites de- 
coupled: 

(cosidi - e i+rCu - $^)) = |(e^)o| 2 (e- 1 *-) 



Since, 



")o = A and finite, we conclude that 
>s(Qi - 9 i+u - 27r$^))o oc h 2 , 



(57) 



(58) 



i. c. finite only in the ordered phase of the dual theory 
(20), i. e. in the disordered phase of the original XY 
model. 
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XV. APPENDIX B 

Here I provide a different derivation of the dynamics of 
the gauge-field a at T ^ starting from the Hamiltonian 
for the Coulomb plasma. Assume a collection of N + (iV_ ) 
vortices (antivortices) at the positions {r,}. The Hamil- 
tonian of the vortex system is 



H v 



1 N 



rj), 



(59) 



where v(r) w — In \ f\, at large distances, and N = N + + 
N~, qi = ±1. The partition function of the vortex system 
Z v can then be written as 



Z v 



oo 



A+!A+! 

N+;~=0 A B A B 



(60) 



fa/2) 
N+\N~\ 



N 



Hdn 



where = A^ (_) + N^~\ and y is the bare vor- 

tex fugacity. The combinatorial factors serve to ensure 
that in Z v one sums over all possible divisions of vortices 
and antivortices into groups A and B, and divides by 
the number of combinations. With this symmetrization 
the symmetry between up and down spin in the original 
Hamiltonian (2) will be preserved in the Dirac theory for 
neutral spinons. This also guarantees that on average 
there is an equal number of vortices (and antivortices) in 
both groups. 

Next, introduce the vorticity densities in Z v by insert- 
ing the unity 

, n a 
1 = / D[ PA }6( PA (f) - J2 HAS(r- r iA )), (61) 
J »=i 



and similarly for B. The gauge field then becomes 
(X7xa(r)) T = n(p A (r)-p B (r)), 



(62) 



in the transverse gauge V • a — 0, and the index denotes 
the t component, v is defined the same way except with 
the plus sign between p A and pb- 

By introducing two auxiliary fields <fr A and $b to en- 
force the constraints, after the integration over the den- 
sities the partition function may be written as 



Z v 



N 



^ ■ n+wx\n+Wb\ 



J D[a, tf, *+,$_] (63) 



exp— [ — \ — f drdr 1 B(r)v(f — f l )B(f l ) 
27HT J 

- ^ In / drexp(iq ia $ a (r))], 



where = $ A ± $b, B(r) = (V x v) T , and b(r) 

V x a) T . Performing the summations yields 

Z v = J D[a, *?,*+,$_] (64) 

exp — [ - ^ J drdr 1 B(r)v(r — r')B{r') 

— y / dr(cos$ A (r) + cos$s(r))]. 



Finally, neglecting the coupling to the charge current, 
the Gaussian integration over v (i. e. B) gives 

Z v = J D[d,<P + ^^cxp-[J rfr[|(V$+) 2 + (65) 

+-b(r)®_(r) - 2ycos($ + (r))cos($_(r))], 

7T 

where I also have rescaled the $ fields by a factor of 
two. The last expression is then analogous to the T = 
expression in the Eq. (18) with x finite and without 
the dual angles B a ^b- By introducing a source term in 
the action, <~ i J j (r)b(r) / V , and integrating over 6, one 
readily finds 



((V x a(f)) r (V x a(f")) T ) = (y)6(f- ?), 



(66) 



where (y) = j/7r 2 (exp(z<i> + )), with the average to be taken 
at $_ = a = 0. One recognizes (y) as the renormalized, 
or running, fugacity in the Kosterlitz-Thouless scaling, 
which signals the appearance of free vortices, (y) plays 
the role analogous to the vortex loop condensate in 2+1 
dimensions, in providing a mass for the field $ + in the 
Eq. (65). This implies the Maxwell term at T ^ for 
the r component of V x a once fluctuating vortices are 
integrated out. 



XVI. APPENDIX C 

For completeness, here I outline the derivation of the 
result that chiral symmetry in isotropic massless QED3 
is spontaneously broken for N < N c , with N c finite, at 
any value of the coupling constant. 

Rescaling the momenta p/m — > p and self-energies 
T,(p)/m — > m and H(p)/m 2 H(p), after taking the 
limit q — » in the Eq. (32) we find 



1 



ir 2 m J a 



p 2 np) 



(p 2 + ^(p)W + Ii{p)Y 



where the polarization is now 

n(p) = ^F^/(p), 



47rm 



(67) 



(68) 



ia,a=A,B 



with 
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f(p) = (2 + 



P 



(69) 



to the leading order in N |20|| . We see that the right-hand 
side of the Eq. (67) is a decreasing function of m, so for 
rn ^ solution to exist we just need RHS(m — 0) > 1. 
This is satisfied for N < N c , where 



dp 



P 2 n P ) 



(p 2 + v 2 (p))f(pY 



(70) 



As defined, S(0) = 1, and one expects S(p) to vanish at 
large momenta. Also, f(p) w irp/2 for p >> 1, so the 
integrand at large argument behaves like ~ T,(p)/p. N c is 
therefore finite, and independent of the coupling constant 
($). Its precise value in the large-N approximation will 
depend only on the function S(p) at N = N c , and can be 
obtained by solving the differential equation equivalent 
to the integral equation (67) |20| (see Appendix E). This 
yields N c = 32 /tt 2 , not far from the results of other more 
elaborate computations that go beyond the leading order 



XVII. APPENDIX D 

Here I discuss a different representation of the quasi- 
particle action, more in line with the previous work [Q. 
This should serve to underline the difference between the 
approximate chiral SU C (2) symmetry, and the exact spin 
rotational 5*0(3), also present in dSC. It is only the lat- 
ter that will appear in the different version of the theory 
considered here and in jf) , while the chiral symmetry will 
remain completely obscured. 

I start again from the same quasiparticle action in the 
Eq. (2), but now introduce the four-component field as 

*i t (2)(^ w «) = (c+(fc,w„),c_(-fc, -w„), (71) 
c f _(k,u; n ), -c+(-fc, -u n )). 

By linearizing the spectrum and by retaining only the 
modes near the four nodes, the continuum theory may 
again be written as 

S[V']= J d 2 r dT^' 1 f [d T + M 1 Vfd x + M 2 v A dy]^[ (72) 

+(l->2,a;<-»y), 

but this time with different form of the matrices Mi and 
M 2 : Mi = —il (g) 173, and M 2 = il ® <J\- Introducing 
7o = 03 (g) 02, for example, the theory becomes 

S[tf'] = J d 2 r J drf i[ 7o 9 T + nv f d x + 72V A d y ]% (73) 

with 7i = (T3 (g) <7i and 72 = 0-3 <£> 03 . It is interesting to 
consider the generators of the global U(2) = U(l)xSU(2) 



symmetry per Dirac component present in this represen- 
tation of the theory. They are I4 = I (8 I, 73 = a\ % I, 
75 = —(72 ® I, and 735 = (73 I, respectively. One may 
recognize the U(l) factor as representing now the contin- 
uous translations, since under a translation c a [k, co) — > 
e lk ' R c a (k,ui), the Dirac field now transforms as 



%(f,r) ->e iK *- R 9' i (f+R,T). 



(74) 



The SU (2) operators, on the other hand, are nothing but 
the spin rotations. In fact, the above U(2) is an exact 
symmetry of the Hamiltonian (2), and is present even if 
all higher order derivatives are retained. 

Including the coupling to vortex loops via massless 
gauge field in the above representation of the problem 
then may spontaneously induce only the d+ip insulator. 
This breaks two of the above generators, which then sim- 
ply rotate the spin-axis. Translational symetry is, on the 
other hand, always preserved in this formulation, and the 
SDW remains invisible. 



XVIII. APPENDIX E 

Here I provide the details behind the numerical solu- 
tion of the Eqs. 50-51. Since we are interested only in 
the qualitative effect of the U-term, it will suffice to as- 
sume that the fermion mass is small, m | ($) | 2 , so that 
one can neglect the p 2 term compared to H(p) in the Eq. 
(51), and take 



n(p) 



mm 



-p, 



(75) 



appropriate for p >> m. This approximation is known 
to lead to even quantitatively good result for the mass 
for N as low as unity, when U = |^0| . Performing the 
angular integrals then gives 



£(«) = x 



Nw 2 q Jo 



A dh kZ(k)(k - (k - q)8(k - q)) 



S 2 (/fc) 



(76) 



Differentiating twice one finds that this i nteg ral equation 
is equivalent to the differential equation |BCJ : 



dq {q dq {q>> Ntt 2 g 2 + S 2 (g)' 
with the boundary condition 

A£'(A) + £(A) = x, 

and with 



U 



X : 



(2^) 2 



dq- 



g 2 s(g) 
2 + s 2 (?)' 



(77) 



(78) 



(79) 



Here I take A = 
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The above equations may now be studied by assum- 
ing q 3> which leads to linear equation that can 
be exactly solved (8^]. This yields, for example, the 
well known transition line in the g — N plane: g c (N) 
(1/4)(1 + y/l- (N c /N) 2 , for N > N c , g c < 1/4 for 
N = N c , with N c = 32/tt 2 . To determine the size of E(0), 
however, one needs to solve the full non-linear equation. 
This may be accomplished, for example, by choosing a 
value for Xt assuming S(A) next, and then iterating back 
to find E(g) for < q < A. The solution is found by 
tuning S(A) to achieve S(0) finite. One then computes 
the value of g = UA/(2tt) 2 from the assumed x an d the 
found £((?). This procedure leads to the Fig. 4. 
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